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Abstract. These notes are devoted to a summary on the mean-field limit of large en- 
sembles of interacting particles with applications in swarming models. We first make a 
summary of the kinetic models derived as continuum versions of second order models for 
5^ ' swarming. We focus on the question of passing from the discrete to the continuum model 

Oh| in the Dobrushin framework. We show how to use related techniques from fluid mechanics 

•^T equations applied to first order models for swarming, also called the aggregation equation. 

We give qualitative bounds on the approximation of initial data by particles to obtain 
the mean-field limit for radial singular (at the origin) potentials up to the Newtonian 
singularity. We also show the propagation of chaos for more restricted set of singular 
potentials. 
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1. Introduction 

In the last years, we have seen the development of a great deal of different models in the 
biology, applied mathematics, and physics literature to describe the collective behavior of 
individuals. Here, individuals may mean animals (insects, fish, birds,...), bacteria, and even 
robots. Most of these models involve the nonlocal character of the interaction as a basic 
modelling pillar, see for instance [201 E31 EOj [78]. In fact, one of largest source of collective 
behavior models comes from control engineering. There, the aim is to produce a suitable 
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control of the movement of small squads of robots in order to perform unmanned vehicle 
operations, for instance [72]. Even, these ideas have been proposed to model crowd motion, 
including more "intelligent" particles deciding their movement based on optimization of 
certain quantities: time to exit from a room or a stadium, for instance |19j . 

Either in social or in biological sciences, these models encounter many interesting features 
such as the spontaneous formation of different pattern behaviors. When we talk about 
patterns, we do not mean static patterns like in the study of crystals but rather dynamic 
patterns leading to the collective motion of the individual ensemble. For instance, two of 
the main collective motion patterns studied in different models are the flock and the milling 
behavior, see [37J ESI EH [281 ES] ■ I n the flock pattern, individuals achieve a consensus on the 
direction or orientation towards some objective, producing as a consequence a particular 
spatial shape showing their preferred comfort structure. This kind of swiftly moving flocks 
have been reported in many species although the most spectacular or bucolic ones are the 
bird flocks, starlings for instance. In the mill pattern, individuals arrange into a kind of 
vortex like motion around some point. This particular moving pattern has been observed 
in fish schools. Hundreds of movies can be easily accessed through internet search showing 
them. 

There are many reasons one can argue, why such a large number of individuals react to 
external stimuli producing these macroscopic patterns without seemingly the presence of a 
leader in the swarm. Hydrodynamic enhancement, predators avoidance, social interactions, 
spawning survival rate, and many others have been proposed to explain this behavior in 
different species, see [7Tj . 

One of the main question in describing this behavior by mathematical models is how 
to include the interaction between individuals. In any case, there is a consensus that the 
modelling starts from particle-like models as in statistical physics. These particle models 
are also called Individual Based Models (IBMs) in the community. They are usually formed 
by a set of differential equations of Newton type (called 2nd order models) or by kinematic 
equations where the inertia terms are neglected (called first order models). Essentially, 
by admitting that the inertia term is negligible, we assume that individuals can adjust 
to the velocity field instantaneously, an approximation valid when their speed is not too 
large. In any case, these first order models were proposed in the literature derived in a 
phenomenological manner [671 ESI EH [TBI [77J [38] . The literature on first and second order 
models for swarming has increased exponentially fast in the last few years. Many of these 
models find also their origin in social sciences, where consensus or opinion formation was 
also described in similar grounds. Another typical ingredient in these models is some kind 
of noise leading to systems of SDEs. In this work, we will not discuss how to incorporate 
noise in these models, we refer to [16] and the references therein. 

Most of these models are based on discrete approaches incorporating certain effects that 
we like to call the "first principles" of swarming. These first principles are based on mod- 
elling the "sociological behavior" of animals with very simple rules such as the social ten- 
dency to produce grouping (attraction/aggregation), the inherent minimal space they need 
to move without problems and feel comfortably inside the group (repulsion/collisional avoid- 
ance) and the mimetic adaptation or synchronization to a group (orientation/alignment). 
Even if these minimal models contain very basic rules, the patterns observed in their simu- 
lation and their complex asymptotic behavior are already very challenging from the mathe- 
matical viewpoint. The 3-zone models including attraction, repulsion, and alignment effects 
are classical in fish modelling [H [52] for instance. Based on them, one can incorporate may 
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other effects to render more realistic the outputs of the simulations and the models, see [8] 
for fish schools or [51] for birds flocks. We also refer to the reader to the recent review [27] 
about the kinetic modelling of swarming. 

To the eyes of a kinetic theorist or a statistical physicist, studying such systems of ODEs 
when the number of individuals get large is doomed to failure. Dynamical system approaches 
are quite useful but they typically have huge problems to describe large systems of particles. 
A classical approach to attack the problem is to pass to a continuous description of the 
system. This means to go from particle descriptions to kinetic descriptions where the 
unknown is the particle density distribution in position-velocity (phase) space for 2nd order 
models or in position space for 1st order models. 

Going from particle to continuum descriptions is one of the most classical problems in 
kinetic theory. It is at the basis of the derivation of the mother and father kinetic equations, 
namely: the Vlasov and the Boltzmann equations. A rigorous derivation of the Boltzmann 
equation from the Newtonian dynamics has only been given for short times (of the order 
of the average time of first collision), see |57| |40| . In that case, interactions between the 
particles are modelled by short-range potentials leading to collision kernels. The question 
of the derivation of the Boltzmann equation from particles with jump processes was also 
raised and solved by [53], and further results are given in the recent important work by [65!. 
The derivation of the Vlasov equation is well understood only for regular or not too singular 
potentials P2J [681 EHl EH]. In fact, a full derivation of the Vlasov- Poisson system in 3D is 
also lacking. The problem of passing to the limit from particle to continuum models like the 
Vlasov equation is called the mean-field limit. This name just comes from the fact that the 
resulting equation is a kind of averaged version of the interaction between the large number 
of individuals. Moreover, the resulting equation gives the typical behavior of one isolated 
individual among all the others since they are assumed to be completely indistinguishable. 

Finally, there are other famous mean-field limit equations, such as the Euler and the 
Navier-Stokes equations for incompressible fluids, see [631 [62]. It nas been extensively used 
for numerical purposes that both equations in the 2D incompressible case can be derived 
from particle approximations, called vortex point approximations. The convergence in the 
viscous case has been rigorously proved for very general initial data [691 [39] . In the non- 
viscous case [73] proves that particle approximations converge towards solutions of the Euler 
equation, but they may not converge to the good solution because of the lack of uniqueness 
in the Euler equation, see [35]. However, in the case where the initial particles are equally 
spaced on a grid to approximate a smooth solution of the Euler equation, the convergence 
was shown in [33]. These vortex methods have been proven to be convergent and estimates of 
the error committed have been obtained in recent works using optimal transport techniques 
|48| but not for the real Euler equation in 2D. 

The aim of this work is to show in detail a particular example of the mean field limit 
in the case of first order models not covered in the previous literature. Nevertheless, we 
will first discuss some of these issues for 2nd order models summarizing results in [221 [16] . 
We will also discuss that the spatial shape of the main patterns: flock and mills, are given 
by stationary solutions of the 1st order models. This gives another reason from a more 
conceptual mathematical viewpoint of reducing to 1st order models. Section [3J will be 
devoted to obtain the mean field limit to the so-called aggregation equation for singular 
potentials recovering some of the models studied in [101 lllj . Here, the idea is to assume 
that we have solutions of the model in better functional spaces due to the singularity 
of the potential, but we have to pay in terms of conditions on the initial distribution of 
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particles (how they are distributed) in such a way that the particle solution converges to 
the continuum solution of the aggregation equation as N — > oo. We will make use of 
similar arguments to [48] to show the mean-field limit for first order swarming models with 
singular potentials up to the Newtonian singularity. In Section HI we study a local existence 
of a unique L p -solution for the aggregation equation. This complements the well-posedness 
theory in [11]. Finally, Section [5] is devoted to show the propagation of chaos property for 
the aggregation equation. This property is very important from the physical relevance of the 
kinetic and aggregation models, since it states that one can derive the mean- field equations 
under quite generic randomly generated initial location of the particles. We are only able to 
show it for a more restricted set of singular potentials with respect to the mean-field limit. 

2. The Dobrushin approach 

2.1. Some Individual Based Models. As we described in the introduction, the mod- 
elling in swarming starts by introducing some particle models, IBMs in the jargon of this 
community, incorporating some of the basic effects: repulsion, attraction, and alignment. 
Let us discuss briefly some of these models, starting with the ones that have recently at- 
tracted more attention due to their simplicity while having a rich mathematical structure 
and pattern formation. One of these models was introduced by the UCLA group in |37] 
and it consists in Newton's like equations where all the effect of repulsion and attraction is 
encoded via a pairwise potential W : M. d — >■ R. A popular choice for the interaction potential 
W is the Morse potential given by 

W(x) = -C A e- {x{/iA + C R e-W /iR , (2.1) 

where Ca,Cr and £a,£r are the strengths and the typical lengths of attraction and re- 
pulsion, respectively. They are chosen for having biologically reasonable potentials with 
C = Cr/Ca > 1 and £r/£a < 1> see ]29] for other nice choices of the interaction potentials 
and a deeper discussion on the issue of biologically relevant interaction potentials. Apart 
from this, the other effect included is the tendency of the particles to travel asymptotically 
at a fixed speed as in |59] . Consequently, a term producing a balance between self-propulsion 
and friction is introduced imposing an asymptotic speed to the particles (if other effects are 
ignored), but it does not influence the orientation vector. The resulting ODE system reads 
as: 

^ = (a-p\v i \ a )v i -j;Y ! VW(\xi-x j \), (z = l,...,iV). 

where a, /3 are nonnegative parameters, determining the asymptotic speed of particles given 
by y/a//3. Here, the potential has been scaled depending on the mass of each particle as 
in |25] and in such a way that the effect of the potential per particle diminishes while the 
energy is of constant order as the number of particles N diverges. This scaling is the so- 
called mean-field scaling, see the introduction of [13J for a nice discussion of the different 
searings in first order models. 

Another popular IBM including only the alignment effect is the so-called [33] model. 
Each individual in the swarnr changes its velocity vector based on the other individuals by 
adjusting/ averaging their relative velocity with all the others. This averaging is weighted 
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in such a way that closer individuals have more influence than further ones. For a system 
with N individuals the Cucker-Smale model reads as 

( dxi 
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for some 7 > 0. 

Associated to the above models, one can formally write the expected Vlasov-like kinetic 
equations as N — > 00, see for instance [25], leading to 

dtf + v • V x f - (VW * p) ■ V v f + div„((a - (3\v\ 2 )vf) = 0, 
where p represents the macroscopic density of /: 



(2.2) 



p(t,x) :-- 



f(t,x,v)dv for i > 0, x£ 



The Cucker-Smale particle model leads to the following kinetic equation: 

df 



Ot 



+ v-V x f = V v -[C[f]f], 



(2.3) 



where £[f](x,v,t) = (H * f) (x,v,t), with H(x,v) = w(x)v and * standing for the con- 
volution in both position and velocity (x and v). We refer to [Ml HSJ [23 [26] for further 
discussion about this model and qualitative properties. 

Moreover, quite general models incorporating the three effects previously discussed with 
additional ingredients, such as vision cones or topological interactions, have been considered 
in [271 [60[ [H EJ E] . In particular in [60] , they consider that the N individuals follow the 
system: 

dxi 



dt 



Vi, 



dv i __ F A , F I 
-dJ- Fi +F ^ 



(2.4) 



where F^ is the self-propulsion generated by the ith-individual, while F( is due to inter- 
action with the others. The interaction with other individuals can be generally modeled 
as: 

N N 

Fl = F^ + F^ = Y J 9±(\x t -x \)^l- + Y,h ± i\v l -v 3 \)^l- 

J=l J J=l J 

Here, g + and h + (g_ and /i_) are chosen when the influence comes from the front (behind), 
i.e., if (xj — x{) ■ v-i > (< 0); choosing g + 7^ g_ and h + ^ h- means that the forces from 
particles in front and those from particles behind are different. The sign of the functions 
g±(r) encodes the short-range repulsion and long-range attraction for particles in front of 
(+) and behind (-) the ith-particle. Similarly, h + > (< 0) implies that the velocity- 
dependent force makes the velocity of particle i get closer to (away from) that of particle 
3- 
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Some of these models, for instance [TJ [3j \M\, include sharp boundaries for the vision 
cone or for the interaction with the nearest neighbors. As we shall see later, these are 
typical situations in which the mean-field limit for general measures will not work. By 
sharp boundaries we mean that the functions involved in the kernels such as w(x), g±, or 
h± are given by characteristic functions on sets depending on the location/velocity of the 
agent. 

2.2. Basic tools in transport distances. In this subsection, we present several defini- 
tions of Wasserstein distances and their properties. 

Definition 2.1. (Wasserstein p-distance) Let pi, p 2 he two Borel probability measures on 

R . . Then the Euclidean Wasserstein distance of order 1 < p < oo between p\ and p2 is 

defined as 

l/p 

d P {pi,P2) :=inf ( / \x - y\ p d^(x,y) 



7 

and, for p = oo (this is the limiting case, as p — )• ooj, 

doo(pi,P2):=inf 7- sup \x - y\ , 

7 \ (x,y)eR d xR d J 

where the infimum runs over all transference plans, i.e., all probability measures 7 on B^xIR^ 
with marginals p\ and p 2 respectively, 



and 



J 



4>(x)d'j(x,y) = / 4>(x)pi(x)dx, 

JR d 

(f)(y)d>y(x,y) = cp{y)p 2 (y)dy, 



for allege C b (R d ). 



We also remind the definition of the push-forward of a measure by a mapping in order 
to give the relation between Wasserstein distances and optimal transportation. 

Definition 2.2. Let p\ be a Borel measure on R and T : M — > 1R be a measurable 
mapping. Then the push-forward of p\ by T is the measure pi defined by 

p 2 {B)= Pl {T-\B)) for BcM d , 

and denoted as p 2 = Tffp\. 

The set of probability measures with bounded moments of order p, denoted by V P (M. ), 
1 < p < 00, is a complete metric space endowed with the p-Wassertein distance d p , see [80] . 
We refer to [41| [6T] for more details in the case of the doo distance. 

Remark 2.1. The definition of p 2 = Tjfp\ is equivalent to 

4>(x)dp 2 (x) = / 4>(T(x)) dpi(x) , 



for all cj) € Cb(M. ). Given a probability measure with bounded p-th moment po, consider two 
measurable mappings X\,X 2 : M. d — > M. d , then the following inequality holds. 

dP(X 1 # Po ,X 2 #p ) < [ \x-y\Pdry(x,y) = [ \X x {x) - X 2 (x)\v d Po {x) . 

JM d xR d JR d 
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Here, we used as transference plan 7 = {X\ x X2)#po in Definition \2.1\ 

2.3. A quick review of the classical Dobrushin result. Under smoothness assump- 
tions on the ingredient functions of the swarming models, one can use adaptations of the 
classical result of |36j to obtain what is called the mean-field limit equation for general par- 
ticle approximations of any initial measure. These arguments are classical in kinetic theory 
and were also introduced in [181 168] . making use of the bounded Lipschitz distance, and 
reviewed in [Ml [79], see also [75j[6l] for the case with noise. The bounded Lipschitz distance 
or dual IU 1,00 -norm is equivalent to the Wasserstein distance d± for compactly supported 
measures. This strategy works as soon as the velocity field defining the characteristics of 
the model is a bounded and globally Lipschitz function whose dependence on the measure 
itself is Lipschitz continuous in the d\ sense. These ideas were improved to allow for locally 
Lipschitz velocity fields for compactly supported initial measures in |22| and for suitable 
decay conditions at infinity and with noise in |16j . With these techniques one can include 
quite general kinetic models for swarming in this well-posedness theory. 

Let us introduce some notation for this section: A = P C (R x R ) denotes the subset of 
"P(R x R ) consisting of measures of compact support in R x R . On the other hand, we 
consider the set of functions B := Lip ioc (R x R ), which in particular are locally Lipschitz 
with respect to (x, v). Br will denote the ball centered at of radius R in R x R. 

Let us consider general operators from measures to vector fields, %[■] : A —> £>, satisfying 
the following hypotheses: for any Rq > and f,g G A such that supp/ U suppg C -Br , 
there exists some ball Br C R rf x R rf and a constant C = C(R, Rq) > 0, such that 

||W[/]-H[&]IU«» (Bjl) <Cdi(/,p), (2.5) 

Up R (H[f]) < C, \\M[f]\\ Loo{BR) < C. (2.6) 

Here, Lip^(-) denotes the Lipschitz constant of a function in Br. 

Given / G C([0,T],V c (B Ro )), and for any initial condition (X°,V°) G R d x R d , the 
following system of ordinary differential equations has a unique locally defined solution 
thanks to conditions (12.51) 



—X = V, X(0) = X° (2.7a) 

fv = H[f(t)](X,V), V(0)=V°. (2.7b) 

We will additionally require that the solutions to that system are global. Of course, this 
is a requirement that has to be checked for every particular model. We prefer to give a 
general condition which reduces the problem of existence and stability to the simpler one 
of existence of the ODEs. Under the above conditions, the existence and uniqueness of 
associated transport equation 

d t f + v-v x f-v v -[n[f}f} = o. (2.8) 

was obtained in [22] to which we refer for full details. In [22], the interactions %[f] = 
(a — {3\v\ 2 )v — VW * p and 1~L[f] = H * f corresponding to (j2.2|) and (I2.3p . respectively, and 

H[f] = F A (x, v) + G(x) * p + H(x, v) * /, 

with Fa, G and H given functions satisfying suitable hypotheses, such that the kinetic 
equation (j2.8|) corresponds to the model (j2.4|) are investigated. 
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Theorem 2.1. Given an operator %{■] : A —> B satisfying Hypotheses (|2.5p and (j2.6|) 
for which the characteristics (I2.7ap - (|2,7bj) are globally well-defined, and /o a measure on 
R d x ~R d with compact support. There exists a solution f on [0, +oo) to equation (12.81) with 
initial condition /q. In addition, 

feC([0,+oo);V c (R d xR d )) (2.9) 

and there is some increasing function R = R(T) such that for all T > 0, 

supp ft C B R{T) CK^xR* 1 /or all t G [0, T]. (2.10) 

TTiis solution is unique among the family of solutions satisfying (]2.9h and fl2. 10f) . Moreover, 
given any other initial data go £ ^(M x Mr) and g its corresponding solution, then there 
exists a strictly increasing function r(t) : [0, oo) — > Mq with r(0) = 1 depending only on Ti 
and the size of the support of /o and go, such that 

di{ft,9t)<r{t)di(f ,go), t > 0. 

The stability theorem 12.11 gives in particular a rigorous derivation of the kinetic equation 
(|2.8p from the large particle limit of the system of ordinary differential equations. This 
is the exact statement of the mean-field limit for general measures as initial data. Let us 
consider the system of ordinary differential equations: 

x t = Vi, i = l,...,N, (2.11a) 

v i = Y,m j H\f N {t)]{x i ,v i ), i = l,...,N. (2.11b) 

where mi, . . . , mjv > and ^i 771 * = 1 an< i / ^ s defined next. Under the conditions of 
Theorem 12.1^ we first notice that if Xj, Vi : [0, T] — > W d , for i = 1, . . . , N, are a solution to 
the system (I2TTT) . then the function f N : [0,T] -)■ P c (M d x M d ) given by 

v 

if : =Xl mi ^i(*)^(*)) ( 2 - 12 ) 

is the solution to (12. 8p with initial condition 

/o r = E m * <y (**(o),«i(o))- ( 2 - 13 ) 

8=1 

In fact, the solution f|2. 12j) is called the empirical measure associated to the system of ODEs 
(|2.1ip . We finally write the full statement of the mean- field limit in the Dobrushin strategy. 



Corollary 2.1. Given /o G V C (M. x K ) and H[f] satisfying the conditions of Theorem 
\2.1\ take a sequence of /g of measures of the form (|2.13p (withrrii, x,(0) and Vi(0) possibly 
varying with N ), in such a way that 

lim d 1 (/ iV ,/ ) = 0. 

N— >oo 

Consider / t i/ie empirical measure associated to the solution of the system (|2.1ip to£/i 
initial conditions Xi(0), Vj(0). Then, 

lim ci 1 (// v ,/ t ) = 0, (2.14) 
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for all t > 0, where f = f(t,x,v) is the unique measure solution to eq. (J2.8P with initial 
data /o- 

This section can be directly applied to the models recently introduced in p] to account 
for vision cones and braking/acceleration of individuals and those in [31 27] to include 
topological (nearest neighbours) interactions once the parameter functions are smoothed 
out to avoid sharp boundaries. 



2.4. First-order models: Aggregation Equation. Summarizing the previous subsec- 
tion, under suitable smoothness of the parameters involved in the swarming models, the 
empirical measures are solutions themselves of the Vlasov-like kinetic equation (12. 8ft . Thus, 
an stability result in d\ with respect to the initial data is enough to conclude the mean-field 
limit. Let us consider one of the particular examples in subsection 2.1, the model intro- 
duced in [37] with the Morse potential (|2.1j) . This potential does not satisfy the smoothness 
assumption in Theorem 12. 11 In principle, one cannot expect to have a mean-field result for 
general measures as initial data and for general approximations by particles. In fact, we 
do not have a well-posedness theory for such initial data in those cases. However, one can 
develop well-posedness theories in better functional spaces, say L 1 n L p (W d x M. d ) for the 
initial data and then impose suitable conditions to the distribution of the approximated 
particles initially to be able to conclude the mean- field limit (J2.14J) . This is the strategy 
that have been followed in [19] for the classical Vlasov equation and in |48j for Euler-like 
equations in fluid mechanics. 

In the next sections, the objective is to show this strategy applied to a simpler swarming 
model than the ones showed above. We will showcase these tools in the case of the so-called 
aggregation equation. Let us assume that we have just particles interacting through the 
pairwise potential W(x). Assuming that the variations of the velocity and speed are much 
smaller than spatial variations, see [66], then one can neglect the inertia term in Newton's 
equation to deduce that 



dX 

- Y^ VW(Xi - Xj) in the N ->• oo limit 



dt 



dp 

— + div (pu) = 

u = — VW * p 



Another reason to study this first order equation is that the stationary states of the first 
order model determine the spatial shape of the flock solutions to the second order models, 
see [29j . Let us remark that one could apply the Dobrushin strategy to the aggregation 
equation for C 2 (M. d ) smooth potential with at most quadratic growth at infinity by following 
the same argument as in Theorem (|2.ip . This argument was detailed in a nice summer school 
notes in 



3. Mean-Field Limit for the Aggregation Equation 

Now, we analyse the mean-field limit of the first order model for swarming introduced 
in the previous section. More precisely, we will study sufficient conditions on the initial 
distribution of particles for the convergence of a particle system towards the aggregation 
equation. This model consists of the continuity equation for the probability density of 
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individuals p(x, t) at position x G M. d and time t > given by: 

dtp + V • (pu) =0, t > 0, x G M ', 

«(t,a;):=-VW*p, t > 0, x G K d , (3.1) 

k p(0,x) :=p°(x), x£R d , 

where u(x, t) is velocity field non-locally computed in terms of the density of individuals. 

As an approximation by particles of the aggregation equation (|3.ip . we consider the 
following ODE system: 

j& ( 3 - 2 ) 

J>Q(0)=X°, i = l,...,N. 

Here, {Xj}^ and {m,j}^ 1 are the positions and weights of i-th particles, respectively. We 
define the associated empirical distribution pjv(£) as 

N N 

i"at(*) = y2 m i s Xi(t), y2 m i= po(x)dx = l, (3.3) 

i=l i=l ^ 

with mj > 0, i = 1, . . . ,N. As long as two particles (or more) do not collide, and if we 
set VW(0) = (arbitrarily if there is a singularity), then p^r satisfies (|3.ip in the sense 
of distributions, i.e., pjv(i) and p(t) satisfy the same equation. In this framework, the 
convergence: 

"p° N — *■ p° weakly-* as measures =^- pjv(t) — ^ p(i) weakly-* as measures 
for small time or for every time?" 

is a natural question. If the answer is yes, we say that the continuity equation (|3.1|) is the 
mean-field limit of the particle approximation f|3.2|) . In other words, we can say that the 
continuum nonlocal equation (|3.ip has been rigorously derived from particle systems. 

Because of the singularity in the interaction force, the natural transport distance to use 
is the one induced by the doo-topology. Remark that this distance also allows to understand 
linearized stability of particle systems around singular steady state measures with a ring 
shape in first order aggregation models, see [3 [56]. Actually, a local perturbation of the 
dynamical system (|3.2p keeping the number of particles fixed is obtained by transporting 
the particle to other locations nearby. One could even allow for splitting of the mass into 
different particles, but all of them located in a local neighborhood of the unperturbed parti- 
cle positions. Certainly, sending a small portion of mass very far away from the location of 
one particle is not a doo-perturbation of the atomic measure but it is a d p small perturbation 
for all 1 < p < oo. These ideas have also recently been used in [7] to study local minimizers 
of the energy functional associated to (13. ip . 

Another issue to cope with is that we are dealing with particle systems whose character- 
istics may lead to collisions in finite time. Therefore, we will be able to obtain meaningful 
results only on intervals in which collisions are avoided (although in some particular cases 
we can allow collisions). 

We next introduce several notations that are used throughout the rest of this work to 
compare the distance between a solution p(t) of the continuum aggregation equation (|3.ip 
and the empirical measure /ijv(i) defined by (|3.3|) associated to a solution {Xj}^ of the 
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particle system (|3.2p . The main two quantities appearing in this comparison are the d^- 
distance between p(t) and pi\f(t), and the minimum inter-particle distance: 

»7(*) :=<to(A«iv(<),p(t)), rjm(t):= min (|^(t) - X,(t)|) , (3.4) 

l<«5=7<iV 

with 77 := 77(0) and 77°, := r/ m (0). Our strategy does not take advantage, as we do not 
know how, of the repulsive or attractive character of the potentials, being the proof equal 
for both cases. 

A theory of well-posedness for measure solutions have been obtained for the aggregation 
equation (|3.1j) allowing collision of particles in finite time in [23} I24j . In these works, 
the potential is assumed to be smooth except at the origin, where the allowed singularity 
cannot be worse that Lipschitz and the potential has to be A-convex, see [23] for details. 
This convexity allows for attractive at the origin potentials, but not repulsive, with local 
behaviors of the form \x\ b with 1 < b < 2. In these works, the essential tools that allow to 
get the mean-field limit for more singular potentials that quadratic are based on gradient 
flows in the Wasserstein distance d<i sense as in [2] . The additional dissipation in the system 
of the natural Liapunov functional given by the total interaction energy is crucial to get the 
mean field limit for general measures for a potential behaving locally at like W{x) ~ |x|, 
for instance for the attractive Morse potential W(x) = 1 — e' x '. 

In this work, we want to allow for more singular potentials at the origin as in |10t 111] , 
and thus we need to work with solutions in better functional spaces. More precisely, we 
will work with solutions of the aggregation equation (|3.ip in L°°(0,T; (L 1 n L p )(M. d )) with 
1 < p < 00 to be determined depending on the singularity of the potential. We will use the 
notation 

ll/ , ll(L 1 nLP)(M d ) := IIpIIi + \\p\\p: IIpII := IIPllL°°(0,T;(L 1 nLP)(M d )) > 

where ||p||p denotes the L p (M. d )-norm of p, 1 < p < 00. 

In order to make sense of solutions to ()3.ip . we need the following assumptions on the 
interaction potential: we first fix W(0) = by definition, even if W is singular at the origin, 
and 

|W(x)|<^p and \D 2 W{x)\ < j-j^, V x € R d \{0} , (3.5) 

for —1 < a < d—1. Note that due to the assumptions on W, we can always find 1 < p < 00 
such that (a + l)p' < d, and thus VW belongs to W^f (R d ). 

Our results also apply with minor modifications for interaction potentials of the form 
W := W\ + W% with W\ satisfying assumptions (|3.5p . and VW2 being a global Lipschitz 
function, or even more general satisfying a one-sided Lipschitz (or convexity) condition 
y ■ D 2 W2(x)y < C\y\ 2 for all y G W*. This last generalization is important because it is 
satisfied if W2 = c\x\ a , (0 < a < 2) with c positive. So that any repulsive-attractive potential 
W, see [U [7] for a definition, such that W(x) ~ — \x\ b /b locally at x near the origin, satisfies 
assumptions (|3.5p locally with a = 1 — b. Therefore, our mean-field limit results apply 
to locally repulsive potentials with exponent range 2 — d < b < a < 2 and without much 
restriction on the attractive part at +00, i.e., a > 0. We will discuss further on localizing 
assumptions (|3.5p at the end of this section. Finally, we cannot apply our techniques to the 
Newtonian singularity [12] being the limiting case of our strategy as it was the case for the 
Euler-like models in fluid mechanics studied in |48] , 

We next summarize the results on the existence and uniqueness of solutions to the aggre- 
gation equation (|3.ip . For the local well-posedness of solutions to equation (|3.ip . we refer 
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to [9[ UO^ fT H 158]. In particular, unique solutions for the system (|3.ip were obtained in |11] 
with second moment bounded initial data. More precisely, Bertozzi et al. |11| Theorem 1.1] 
showed that if VW G W l > p ' (R d ) and p° G L p (R d ) n V 2 {R d ), then there exists T* > and 
a unique nonnegative solution to (|3.ip satisfying 

p g c([o, r*], (L 1 n L p )(M d )) n c x ([o, T*],w~ 1 ' p (R d )). 

Unfortunately, one can not directly apply those results for potentials satisfying assumptions 
(|3.5p . We will compliment the results in [11] to show the local existence of a unique solution 
to the system (|3.ip with the interaction potential function W satisfying (|3.5p in Section SJ 
We prefer to postpone the well-posedness theory in order to emphasize the mean-field limit 
result contained in the following theorem, whose proof follows the strategy in |48| . 

Theorem 3.1. Suppose the kernel W satisfies (|3.5[) . and let p be a solution to the system 
(I3~TT) uptotimeT>0, such that p G L°°(0,T; (L 1 nLP)(R d ))nC([0,T],V 1 (R d )), with initial 
data p° G (Vi n L p )(R d ), < a < — 1 + d/p' , and 1 < p < oo. Furthermore, we assume p? N 
converges to p° for the distance doo as the number of particles N goes to infinity, i.e., 

dooipx^ ) -> as N -»■ oo, 

and that the initial quantities rp , 77^ satisfy 

(r,0)d/p' 
lim y,;, , = 0. (3.6) 



Then, for N large enough the particle system (13. 2\i is well-defined up to time T , in the sense 
that there is no collision between particles before that time, and moreover 

pN(t) — *■ p(t) weakly-* as measures as N — > 00, for all t G [0, T]. 

Remark 3.1. Let us first discuss the assumptions on the initial data in Theorem \3.1[ The 
mean-field limit is valid for particular approximations p° N of p°, that is, for well chosen 
particle approximations of the initial data. In fact, a procedure to construct initial atomic 
measures approximating the initial condition in the sense of (|3.6p is the following: define 
a regular mesh of size e and approximate p° by a sum of Dirac masses p° N located at the 
center of the cells such that the mass at each particle is exactly equals to the mass of p° 
contained in the associated cell. In that case, we have rp ~ e and rj^ ~ e (for the last 
condition we need that the mesh has some regularity). In that case, the assumption (I3.6P 
is automatically fulfilled since (1 + a)p' < d. Notice that no bound on the masses mi of the 
particles is required. 

Proof of Theorem \3.1\ The proof of Theorem 13.11 is divided into three steps: 

• In Step A, we estimate the growth of the doo Wasserstein distance between the 
continuum and the discrete solutions r] that involves rj itself and rj m in the form: 

^CrjWpWll + rjVP'r,-^). (3.7) 



dt ~ ~" ir " V" ' ' ,m 
In Step B, we estimate the decay of the minimum inter-particle distance rj m , which 
also involves the terms 77 and r\ m in the form: 

^i>-C^||p||(l + ^/^ 1+ «)). (3.8) 

In Step C, under the assumption of the initial approximation (|3.6p . we combine ([37 
and (13.81) to conclude the desired result. 
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Step A.- We first introduce the flows generated by the two velocity fields: u(x, t) = 
— VW* p and u^ := — VW* p^. Let us remark that the convolution in the definition of u^ 
is just a notation for the right-hand side of (J3.2J) since the convolution of a Dirac Delta with 
a (possibly) singular potential is not well-defined. These flows ^n, ^ '■ ^+ x ^+ x K — > K 
are defined as solutions of 

— (*(t; s, x)) = u{t; s, *(t; s, x)), (3 g) 



\l/(s; s,x) = x, 



for all s,t E [0,T], and 



— (^ N (t; s, x)) = u N (t; s, *jv(t; s, a;)), , g ^ 

*at(s;s,3;) = x, 

for all s,t £ [0,7^]. Notice that the solution Xj(£) to the system (|3.2|) is well-defined 
and continuous by the Cauchy-Lipschitz theorem as long as there is no collision between 
particles. Since rf^ > 0, there exists T N > such that rf m (t) > for t G [0, Tq ] by 
continuity. Then the flow map \Pjv(£; s, x) solution to (|3. 10[) is well-defined for t,s G [0, T^]. 
Now, let us check that the flow for the solution associated to the continuum equation in 
(|3.9p is well-defined. Assumptions (J3.5J) imply that 

\VW(x)-VW(y)\< n~, zi ■ (3.11) 

1 v ; vy;i - min(|x|,|y|) a+1 v ; 

One can see this by integrating along a straight line joining x and y but avoiding the 
singularity using a small circle if needed, see [38]. The estimate (|3.1ip implies that the 
velocity field is Lipschitz continuous with respect to the spatial variable. Actually, one can 
estimate it as 



\u(t,x)-u(t,y)\< / \VW(x-z) - VW(y- z)\p(t,z)dz 
Jm. d 

<2\x-y\ ——- n -——p(t,z)dz 

J Rd wm(\x - z\,\y - z\) a + x 

< 4|x - y\ sup / : —- T p(t,z)dz. 

xeR d JR d \ x ~ z \ 

Now, splitting the last integral into the near- and far-field sets A := {z : \x — z\ > 1} and 
B := M. d — A and estimating the two terms, we deduce 

-p(t, z) dz < Hp^IIx + ( f J— -- dy 



\x-z\ a+lHK ' ' - ii^wiu \] B \x-y\^°)v' 

<C\\p\\, (3.12) 

for all x £ M. d due to the assumption (l+a)p' < d. Putting together previous inequalities, we 
get the desired Lipschitz continuity of the velocity field with respect to x, which is moreover 
uniform in time. A similar estimate using (J3.5H shows that the velocity field is bounded, and 
then the flow ^ in (|3.9p is well-defined. Our first aim is to find an expansion of the velocity of 
the doc Wasserstein distance. The idea is similar to the evolution of the euclidean Wassertein 
distance in [301 ED ED] • Fixed < to < min(T, T ) and choose an optimal transport map 
for doo denoted by T° between p(to) and //jv(*o); A*Ar(*o) = 7~°#p(£o)- ^ ^ s known that such 
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an optimal transport map exists when p(to) is absolutely continuous with respect to the 
Lebesgue measure [32]. Then it follows from Theorem 14.11 that p(t) = ^(t;to, ■ )#p(to) and 
obviously /ij\r(i) = ^N(t',to, ' )#A*iv(*o) f° r t > *o- We also notice that for t > to 

T t #p{t)=p N {t), where T* = V N {t\t Q , ■) o T° o V(t ;t, ■). 

By Definition 12.11 of the d p Wasserstein distance, we get 

dP p {p N {t),p{t))< J \*(t;t ,x)-y N (t;t ,T°(x))\ p p(t ,x)dx. 

In the case of p = oo, we obtain 

v(t) = d 00 {p N {t),p{t))< ||*(t;to,0 -* N {t;to,-) oT°||oc. 
We notice that 

[9 N {t',to,7*(x))-V(t;to,x)) = u N (t ,T°(x)) - u{t ,x). 



dt 



t=t 



Thus, writing the integral form, dividing by t — to, and taking the limit t — > t$ we easily 
get 

— \\V N (t; t , •) o T° - *(t; to, -)lloo| + < ||t*iv(*o, >T°- u{to, -)||oo- (3.13) 

We now note that 

UN(t ,T°(x)) -u(t ,x) 

VW(T°(x)-y)dp N (to,y)+ I VW(x - y)p{t ,y)dy 



{VW(T°(x) - T°(y)) - VW(x - y)) p(t ,y)dy. 

For notational simplicity, we omit the time dependency on to in the next few computations. 
This yields that f|3. 13[) can be rewritten as 

d + rj 



dt 



<Csup / \VW{T{x)-T{y))-VW{x-y)\p{y)dy. (3.14) 

x£R d JR d 



x 

Bit 



We decompose the integral on R into the near- and the far-field parts as A := {z : \x — z\ > 
4f]} and B := R d - A as 



B (3.15) 

Z 1 +Z 2 . 



[ \VW(T(x) - T(y)) - VW(x - y)\ P {y)dy =/•••+/ 

JR d J A JB 

For the estimate in the set A, we use 

\T{x) - T{y)\ > \x - y\ - \T(x) - x\ - \T(y) - y\ > \x - y\ - 2rj > — - — 

together with (|3.1ip and (13. 121) to obtain 

T . f 2(\x-T(x)\ + \y-T(y)\) 

11 ~ J A min(|* " y|, |T(x) - T(y)\)^ P{mV 

~ 4r] L (j^vW + J^vw) Piy)dy ~ Cr] L W^yW p{y)dy 
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<C V [ l —- T p{y)dy < Crj\\p\\. (3.16) 

jRd \x - y\ a+L 

For the second part I2, we estimate separately each term using (j3.5[) to deduce 

X 2 < [ / {V \ dy + [ P Mdy 
JB \ x ~ V\ JB Vm 

* {Lw^ dy T' M ' + vM uy T »* <3,17> 

< C (,fl''-" + ifVihTsMr < Ci^'"'- + n d/ ''r,-°)]\p]\ ■ 

Notice that \T{x) — T(y)\ > f] m by definition of the minimum inter-particle distance (|3,4p 
as soon as T(x) 7^ T(y), VW(T(x) — T(y)) = otherwise. 

Finally, we choose two indices i,j so that \Xi — Xj\ = r/ m , then we observe that the 
middle point between Xi and Xj has to be transported by T to either Xi or Xj, and thus 
Vm < 277. Hence by combining (|3.14p - (|3.17p and being to arbitrary in [0, min(T, T^)), we 
have 

- C V \\p\\ (l + v dlv '~\m) < C V \\p\\ (l + ^- {1+q) ) , (3.18) 



< 



di 
for all f G [0,min(T,T^)). 

Step B.- We now focus on showing the lower bound estimate of r] m to make the system 
(I3.18P closed. We again choose two indices i, j so that \Xi — Xj\ = r/ m . Neglecting the time 
dependency to simplify the notation, we get 



— \Xi -Xj\> -\u N {Xi) - u N (Xj)\ 

> - f \VW(Xi -y)- VW(Xj - y)\ dp N {y) 

JR d 

= - [ \VW(Xi - T{y)) - VW(Xj - T{y))\ P {y)dy , 

JR d 

where T is the optimal map satisfying /U/v(i) = T#p(t), for each t G [0,min(T, T^)). 
Similar to (|3.15p . we split in near- and far-field parts the domain R as A := {y : \Xi — y\ > 
277 and \Xj — y\ > 2r/} and B :=M. d — A. We can again use (|3.1ip to deduce 

/ \VW(X t - T{y)) - VW(Xj - T(y))\ p(y)dy (3.19) 

J A 

~ J A ^H\X l -T{y)\,\X--T{y)\)^ p{y)dy 

< 2^\ Xl - X 3 \ ! ( l +1 + 1 +1 ) P (y)dy < C Vm \\p\\, 



where we used that \Xi — T(y)\ > \Xi — y\ — i] > ^\X{ — y\ and similarly for Xj together 
with (I3.12p . For the integral over B, we use that as soon as X, 7^ T(y), then we obtain from 
(|33D that 
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and VW{Xi — T(y)) = otherwise, and similarly for Xj. A simple Holder computation as 
in (|3.12j) implies that 

I p(v)dy<CrfV\\p\\, 

JB 
from which we infer that 

/ \VW(X t - T{y)) - VW{X j - T{y))\ p{y)dy < Crj^' ^\\p\\. (3.20) 

JB 

Putting together (13.191) and (13.201) . we finally conclude that 

n >-Cv m \\p\\(l + V d/p '^ 1+a) ), (3.21) 



dt 
for all t€ [0,min(T,T(f)). 

Step C- Until now, we have proved from (J3.18P and (|3.21|) that 

d + r] 



dt 



<Cr,\\p\\(l + V d /P' V ^ 1+a) ), 

^h. > -Or, II nil (\ + n d /p' n~ {1+a) 
, il ^'ImWPW I r i^ '/ '/m 



(3.22) 



for t G [0,min(T,T Ar )). We first notice from (|3T22]> that if r/ d / p ' r] m {1+a) < 1, then 

r/(i) < 7]°e 2||p||t and r, m {t) > r? m e~ 2Mt t e [0,min(T,I^)). (3.23) 

We now show that (|3.23p holds for time t £ [0, T] when iV goes to infinity, in other words 
that T < T^ when N is sufficiently large. For this, we set 

/(*):=#, g(t):=VZ& and fr := foWi&)-< 1+ «>. 

Note that £at depends on the number of particles N as in (|3.4h . It yields 



<c||p||/(i + ^/ d / p V (1+a) ), 

>-C||p|| 5 (l + Z N fWg-( 1+a '>). 



dt 
dg 

dt 

Since /(0) = g(0) = 1 and £tv — > as N goes to infinity, we obtain that there exists a 
positive constant T^ '(< T^) such that 

t N f d/l/ 9- {1+a) < 1 for te[0,7f], 
for sufficiently large TV. Then it follows from (|3.23p that 

f(t) < e 2Mt and g(t) > e ' 2Mt . 
This yields £ N f d/p ' 9~ {l+a) < ^ N e 2( - d / p ' +( - 1+a ^ p ^, that is, 



iNf d,p 'g~ {l+a) < 1 holds for t < 



MZn) 



2{d/p' + {i + a ))\\ P \y 

so that 



2(d/p' + (l + a))||p|| 
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On the other hand, our assumption for the initial data (|3.6p implies 

v • f T JV N r ln(6v) 

hm ml Tl > hm . - , , ; .... .. = oo , 

TV^oo - Af^oo 2{d/p' + {l + a))\\p\\ 

and thus for N large enough, T < T + < Tq . This completes the proof. □ 

Remark 3.2. One can use almost the same argument with the above to obtain an stability 
estimate in doc,: let p± and p2 be solutions given by Theorem ^. l\ to the system (13. ip satisfying 
(|3.5p . then we have 

-^d 0O (p 1 (t),p 2 (t)) < Cmax(||pi||, ||p 2 ||) d 00 {pi{t),p 2 {t)) . 

In fact, the estimate of mean field limit in Theorem 13.11 holds for — 1 < a < without 
any condition on rp and 77^. This is coherent with the results in |23j in which the mean 
field limit is obtained for all measure initial data without restriction in the way initial data 
are approximated by Dirac masses at least for attractive potentials. 



Corollary 3.1. Suppose the interaction potential W satisfies (|3.5p with — 1 < a < 0, and let 
p be a solution to the system d37TJ) such that p G L°°(0,T; (L 1 nL p )(R d ))nC([0,T],Vi(R d )). 
Suppose that 

doo(VNiP°) ~ y as N — y °°- 
Then for any solution of the ODE system (|3.2|) the associated empirical distributions p,N(t) 
converge toward p(t) uniformly in time: 

sup doa^N^), p(t)) — > as N — > 00. 
te[o,T] 

Remark 3.3. It is remarkable that even if we do not have uniqueness of solution of (13.2|) 
under assumption ()3.5[) with — 1 < a < 0, we get the mean field limit without restriction. If 
one collision occurs, then uniqueness may lost, but the existence of solution is still guaran- 
teed. Thus Corollaru \3.1\ is interesting because it is valid for density solutions to (|3.ip even 
if collisions occur and uniqueness is lost at the particle level. 

Proof of Corollary \3.1l We first notice that the existence of solutions to the ODE system 
(|3.2p is guaranteed thanks to Cauchy-Peano-Arzela theorem since a is strictly negative with 
(13. 5p implies that VW is continuous. One can use the same arguments as in the proof of 
Theorem 13.11 to find 

^ < C sup ( I + [) \VW(T(x) - T{y)) - VW(x - y)\ p{y)dy 

:= /Ci + /C2, 

where the same notation for the sets A and B is used and the time dependency has been 
avoided for simplicity. Using (I3.16P we estimate K.\ by Cr/||/9||. To estimate /C2, we use that 
a < to get 

\VW(T(x) - T(y)) - VW(x - y )\<<L + _^_ , 

rj a \x — y\ a 

and to obtain by Holder's inequality that 

IC 2 < C f ^-dy + £ / p{y)dy < C^'- Q ||p|| p + Crfl^n^ \\p\\ p 
J B \x-y\ a rj a J B 

< Crj d / p '- a \\p\\ . 
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Hence, we have 

^<c v \\ P \\{i + v ^'-^ 

and this yields for sufficiently large N 

r){t) < f^^-(d/p'-a) e -c\\ P \\(d/p'-a-i)t + e ~c\\ P \\(d/ P '~ a ~i)t _ A'd/p'-a-i ^ 

for all t G [0, T]. Note that d/p' — a — 1 > and then, the right hand side of previous 
estimate goes to zero as N goes to infinity This completes the proof. □ 

We next show that there is no collision between particles when the initial quantities rp 
and rftn in (13, 4|) satisfy 

( 0\d/p'-a 

lim ^-L. = 0. (3.24) 

Note that the same strategy as in Remark 13.11 allows us to find suitable approximations for 
the initial data satisfying (|3,24p . 

Corollary 3.2. Under the assumptions of Corollary \3.1\ with — l<a<0,ifwe further 
assume that r/°, n^ satisfy (|3.24p . Then we have that for N large enough, the particle 
system (J3.2J) is uniquely well-defined till time T in the sense that there is no collision between 
particles before that time, and the convergence 

sup d 00 (//7v(t),p(t)) — > as N — )• 00 , 

*G[0,T] 

holds. 

Proof. The proof of Corollary 13.11 shows that for sufficiently large N 

V < (( T ] ) 1 ~ Wp '~ a) e' CMWp '- a - 1)t + e ~C\\p\\(d/p'-a-l)t _ 1 \-d/,'-«-i _ 
For the estimate of rj m , one can obtain from the proof of Theorem 13.11 that 

^ > -C Vm \\p\\ (l + V d/p '- a ^) for all t G [0, min(T, T»)), 

where T^ denotes the first collision time between particles. Then we conclude the desired 
result employing the same arguments in Step C of Theorem 13.11 using f|3.24j) . □ 

As a corollary of Theorem 13.11 we consider interaction potentials under weaker assump- 
tions than ()3.5|) : there exists R > such that W satisfies 

\VW(x)\<^-, and \D 2 W(x)\ < -^, VieB(0,fi), (3.25) 

where B(0,R) := {x G M. d : \x\ < R}. Then one can assume that the initial data p° has 
compact support, and show that the local solution p(t) has compact support on a small 
time interval [0, T\. This is possible since characteristics are locally in time well defined and 
the velocity is uniformly bounded under the assumptions (J3.25P initially. This argument 
was made rigorous under stricter assumptions on the local behaviour of the interaction 
potential but allowing growth of the potential at infinity in [5]. Thus, one can cut-off the 
potential outside a large ball in such a way that the solution is unaffected but the potential 
satisfies the global assumption VW G W 1 * (R d ) entering the well-posedness theory in |llj 
or satisfying (|3.5[) allowing for the application of Theorem 14.11 Concerning the interaction 
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potential W satisfying (|3.25p . the same results of convergence in Theorem 13. II and Corollary 
can be obtained. We leave the details to the reader. 

4. Local existence and uniqueness of L p -solutions 

In this section, we provide a local existence and uniqueness result of weak solutions in 
L p -spaces to the system (|3.ip under the assumptions (|3.5|) . 

As we mentioned before, we can not directly apply the arguments in [11] for the potentials 
satisfying (I3.5p . Of course, we can overcome these difficulties using the property of compact 
supports on the initial data p° (see the paragraph below Corollary 13. 2p . However, we use 
the arguments of dividing near- and far-field parts of the interaction potential function W to 
establish the local existence of a unique L p -solution to the continuity aggregation equation 



Theorem 4.1. Assume that W satisfies the condition (|3.5p . for some < a < — — 1, and 

that p° E Vi (R d ) n LP(M. ), 1 < p < oo. Then there exists a time T > 0, depending only 
on ||p°||p and a, and a unique nonnegative solution to (|3.ip satisfying p E L oo (0,T;L 1 n 
U>(R d )) n C([0,T],Vi(R d )). Furthermore, the solution satisfies that there exists C > 
depending only on \\p \\ p and a such that 

\\p(t)\\ P <C for all t E [0,71. (4.1) 

The velocity field generated by p, given by u = — VW*p, is bounded and Lipschitz continuous 
in space uniformly on [0, T] , and p is determined as the push-forward of the initial density 
through the flow map generated by u. 

Moreover, if pi, i = 1,2, are two such solutions to (13. ip with initial conditions p® E 
ViiW 1 ) n L p (M. d ), 1 < p < oo, we have the following stability estimate: 

—di(t)<Cmax(||pi||, H/02 ||)di(t), 
where d\(t) := di(pi(t) , p 2 (t)) . 

Proof. Let us start by proving the uniqueness. Given two weak solutions pi £ L°°(0, T; L 1 n 
L p (R d ))nC([0,T],T J i(R d )), i = 1,2, to the continuous aggregation equations (|53]1 . consider 
the two flow maps \&j : M. + x R + x M. d — > M. d , i = 1, 2, generated by the two velocity fields, 
i.e., 

— (*j(i; s, x)) = Ui(t; s, *i(t; s, x)) , 
at 

^i(s;s,x) = x, 

where Ui := — VW * pi, t,s E [0, T] and x E Mr. We know that the solutions are constructed 
by transporting the initial measures through the velocity fields pi = ^j#p^, i = 1, 2. 

Let T° be the optimal transportation between /?i(0) and p2(0) for the (ii-distance. Then 
we define a transport (not necessarly optimal) between pi(t) and P2(t) by 

T*(x) = * 2 (t; 0, x) o T°(x) o ^i(0; t, x), r*#A>i(t) = p 2 (t), 

and ^d\(t) < Q(t), where di(i) := di(pi(t), p2(t)) and 

Q(t) := / | W(T*(:r) - T*^)) - VW(x - y)\ Pl (t, x) Pl (t, y)dxdy , 
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where we have used a similar argument as in Step A of the proof of Theorem 13.11 To 
simplify the notation, let us not make explicit the dependence on time. Note by symmetry 
that 

:= Ji + J 2 - 

Straightforward computation using the near- and far-field decomposition as in (|3.12j) shows 
that 

Ji = *[ \T(x) - x\ Pl (x) ( f P2{y) , 1+a dy)dx 

jRd \Jm d \T(x) -y\ L+a ) 

<C\\p 2 \\ [ \T{x)-x\ Pl {x)dx = C\\ P2 \\d 1 {t). 

Similarly using again ([3.12|) . we have Ji < C||pi||di(t). It yields that 

— di(t)<Cmax(||pi||, ||p 2 ||) di(t), 
from which we conclude the uniqueness part of the statement. 

Let us now show the existence of weak solution. Let e > and 9 be a standard mollifier: 

9 > 0, 6 e Cg°(M d ), supp 6 C B(0, l), / 0(a;)da; = 1, 
and we set a sequence of smooth mollifiers: 

We first regularize VTV such as VW £ := (VW) *^ e . Then since VW e is a globally Lipschitz, 
we can apply the theory of [18, 36, 58] which says that there exists a unique global solution 
p e to the following system 

d t p £ + V • {p e u £ ) = 0, t > 0, x£l d , 

Ueftx) :=-VW £ *fe t>0, a;GlR d , (4.2) 

{ p e (0,x):=p°(x), x£R d , 

A standard calculation, see [11] . implies that 

-^WPeW^nLP < C\\p e \\ L i nLP , (4.3) 

where C is an uniform constant in e. Note that the inequality (|4.3p holds only formally for 
the non regularized problem, but it is fully rigorous for the regularized one with W £ . This 
yields that the time of blow-up depends only on the initial data, more precisely ||p°||, and 
not on e. Thus, there exists a T > such that 

sup ||p e || < oo. (4.4) 

£>0 

It follows from (|4.4|) and the evolution in time of the first momentum of p, that this first 
moment is also uniformly bounded: 

sup||xp £ || L oo( 0iT;L i( R d)) < C, 

£>0 
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where C depends only on T, ||x/3°||i, and ||p°||. We leave the details to the reader. Next, 
we show an estimate on the growth of the d\ distance r) £>£ i(t) := di(p e (t),p e >(t)) between p £ 
and p £ >, for e,e' > 0: 

j f r]e,e'(t) < Cmax(||p £ ||, ||/v||) (Ve,e'(t) +e + e'), (4.5) 

where C is an uniform constant in e and e' . We remark that the above estimate (|4.5[) implies 
that {/> e } £ >o is a Cauchy sequence in C([0,T],Vi(M. d )). 

Let us remark that the weak solutions to the regularized problems (|4.2p can be written in 
terms of characteristics. This is a consequence of the fact that the associated velocity field u t 
is bounded and Lipschitz in space, unifromly in time and some standard duality arguments. 
This strategy is explained in detail at the end of the proof of the present Theorem applied 
to the solution of the original problem, and we refer the reader there for details. Since 
solutions are constructed by characteristics, for the proof of (|4.5p we can proceed as in the 
part of uniqueness. Therefore, not making explicit the time dependency, we get 

T+ r,e,e'(t) < / |W E (7» - T(y)) - VW £/ (x - y)\ p £ ,(x)p £ ,(y)dxdy 

< I \VW £ (T(x)-T(y))-VW £ (x-y)\p £ ,(x)p £ ,(y)dxdy 



+ / \VW £ {x-y)-VW £l {x-y)\p £ ,{x)p £ <{y)dxdy 

■/R d xR d 
:=/Ci+/C 2 , (4.6) 

where T is the optimal transportation between p e i(t) and p e (t) for the di-distance. To 
estimate /Ci, we notice that 

(y) ^ f Uy) 



\vwe(x)\< i , n , Uv L a dy+ L 

J\y:\y\<^\\ x -y\ + J\v-\y 



<{y:\y\<^} \*-y\ 1+a " ' J{y.\y\> l 4} \ x ~ V^^ 
ol+« r r Q („,\ 

- w^ I eA ' J)dy + 1(l * 2! » 4 »>,„„ w^w** 

< _C_ C^f 6,(y) C 

— |™|1+Q |™|l+ar L ^ , ,, |t_7/|1+o y ~ I -t-II+Q: ^ ' 

l x l l-^l J {y: e>\y\} l x i/l \ mij \ 

Then we now use again the decomposition (|3.12p as in the part of uniqueness to find 

/C^CmaxtllftlUI/vlDr/^), (4.8) 

where C, \\p £ \\, and \\p e '\\ are uniformly bounded in e and e' thanks to the estimate (|4,4D . 
For the estimate of /C2 1 we claim that 

Cr 
\V(W-W £ )(x)\<-^, (4.9) 

where C is independent on e. 
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Proof of Claim: It is a straightforward to obtain 

\VW £ (x)-VW(x)\< I \VW(x-y)-VW(x)\6 £ (y)dy 

Jj ( i +_• W toM „ («») 

:=£l + £ 2 . 

Noticing that the mollifier properties allow to gain an e factor in front of the integrals, we 
can estimate Ci , i = 1,2 as follows 

A < r^7 / e e (y)dy 



™il+o / , c^/ » |™|l+a' 

X\ ./rod X 



£2<2£ / «^*< * , <4 ' U) 

jRd |x — y| \x\ 

where we used a similar argument to (14, 7p for £2- We now combine (I4,10p and (I4,lip to 
have the inequality (14, 9p . Then we use (14, 9p together with (|3,12p to find the estimate of IC2 

K 2 <C{e + e') f P£ f X)P ^ V) dxdy < C{s + e>) \\p £ , \\ . (4.12) 

jR d xR d \x-y\ 1+a 

This completes the proof of the inequality (|4.5|) by combining (|4.5p , (|4.6p , (|4.8p , and (|4.12p . 
Since p £ is a Cauchy sequence in C([0,T],V\(M. d )), it converges toward a limit curve of 
measures p G C([0, T], Pi(M d )), and we also have p G L°°(0, T; L l r\LP{R d )) from the uniform 
bounded estimate (|4,4p . It remains to show that /9 is a solution of the aggregation equations 
(ETTD . Choose a test function 0(i,x) G C~([0,T] x R d ), then p e satisfies 

p°(x)cj) (x)dx= p £ (T,x)<j)(T,x)dx+ / p £ {t,x)d t <j){t,x)dxdt (4.13) 

JM d JR d Jo JR d 

-f f I p £ (t,x)p £ (t,y)VW £ (x - y) -V<j)(t,x)dxdydt. 

JO JR d JR d 
The first two terms in the rhs of f|4. 13|) converges to 

p(T,x)(j)(T,x)dx + / / p(t,x)dt(f){t,x)dxdt, 
Jo JR d 

since p £ — > p in C([0,T],Vi(R. d )). For the third term in the rhs of (|4.13p . we use the 
estimates P~9|) and (fO|) to find 



II 



p £ (t, x)p £ (t, y) (VW £ (x -y)- VW(x - y)) ■ V(f>(t, x)dxdydt 
'o. 

as e — > 0. It remains to show that 

/ / p £ (t,x)p £ (t,y)VW{x -y)-V<p(t,x)dxdydt 

'0 JR d JR d 

/ / p(t,x)p(t,y)X7W(x-y)-X7(f>(t,x)dxdydt, 
'0 JR d JR d 

as e — > 0. For this, we introduce a cut-off function xs G C£°(]R) such that 

, s _ / 1 if |x| < <5 
XtfW ~ 1 if Id > 25 ' 



0. 
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Then it follows from the weak convergence that 

Pe(t, x)p £ (t, y)(l - xs(x ~ y))VW{x - y) • V0(i, x)dxdydt 

d 

T # 

p(t, x)p(t, y)(l - xaO - y))VW(x - y) • V</>(£, x)dxdydt, 

as e — )• 0, since (1 — x<5( x ~~ y))VW(x — y) ■ V(p(t,x) is a Lipschitz function. We estimate 
the remainder as follows: 

/ / Ps(t,x)p e (t,y)xs(x - y))VW(x - y) -V<j)(t,x)dxdydt 

JR d JR d 

T r ! 



<CS / -—p £ (t,x)p E (t,y)dxdydt 

JO J {(x,y)£R d xR d : \x-y\<2S} \ x ~ V\ 

<CT5\\p £ \\^0 as <5^0. 



Similarly, we have 

rT 



lim 



p(t, x)p(t, y)xs(x - y)VW(x - y) ■ V(f>(t, x)dxdydt 
o 



0. 



Hence, we conclude that p satisfies 

p°(x)(j) (x)dx=-- I p(T,x)(j)(T,x)dx+ p(t,x)dt<j>(t,x)dxdt (4.14) 

R d JO JR d 

- I I I p(t, x)p(t, y)VW(x - y) ■ V0(t, x)dxdydt, 
JO JR d JR d 

forall<£GC£°([0,T] x R d ). 

Now, We notice that a weak solution in p G L oo (0,T;L 1 n L p (M d )) to fl3J]) under the 
assumptions (|3.5p has a well defined flow by using the same arguments as the ones at the 
beginning of Theorem 13.11 In fact, the velocity field is bounded and Lipschitz continuous 
in space with 

\u(t,x) -u(t,y)\ < C\\p\\\x-y\ 
for all x, y G M. d and t G [0, T]. Thus, the flow map 

— (V(t;s,x)) =u{t;s,-$(t;s,x)), 

dt 

^(s; s,x) = x, 

for all s, £ G [0, T] is well-defined. Choosing as test function in (|4.14|) (f)(t, x) = ip(^(t; T, x)) 
for any T G (0, T] with ip G C£°(M d ), it is a straightforward to check, due to the definition 
of the flow map, that 



p°(x)ip(V(0;T,x))dx= / p(T,x)ip(x)dx, 

jR d 

for all ip G C£°(]R rf ), and thus by a density argument we conclude p(T) = ^(T;0,- )#p°. 
Since this argument can be done for all < T < T, this completes the proof. D 
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5. Propagation of chaos 

In most practical purposes to approximate the continuum model by particle systems, it 
is naturally expected that initial positions and velocities will randomly and independently 
be selected. We will show that the empirical measure at time is then close to p° with 
large probability in suitable weak norm. 

In a seminal article [53] , the propagation of chaos was introduced by Kac giving a proof 
for a simplified collision evolution process. He showed how the limit of many particles 
rigorously follows from the property of propagation of chaos. For a classical introduction 
to these topics, we refer to [75]. Later, this property has been studied and developed in 
kinetic theory, [Ml E3 El EQJ [65] . 

Let us introduce the notion of propagation of chaos. Let us consider p N (t,x\, ■ ■ ■ , xjv) 
being the image by the dynamics to the coupled system (13. 2D with iV-equal masses particles 
of the initial law (p )® . We define the /c-marginals as follows. 

Pk(t,xi,--- ,x k ) := / p N (t,x)dx k +i'" ,dx N . 

J R d(N-k) 

Let us choose the initial positions X ,0 := {Xj > }fL 1 as independent identically distributed 
random variables (in short iid) with law p. We can construct the associated empirical 
measure as in (|3.3p by 

1 N 

1=1 

but now understood as a random variable with values in the space of probability measures. 
The propagation of chaos property is defined as follows: for any fixed k £ N, 

Pk ~^ (p)® weakly-* as measures as N — > oo. 

It is classically known |75] that it is sufficient to check this property for k = 2 to derive 
the propagation of chaos. In fact, this is based on the fact that propagation of chaos is 
equivalent to show that the empirical measures /Uj\r(i) converge in law towards the constant 
random variable p(t). 

Theorem 15.11 gives a quantified version of the convergence in probability of pat(£) towards 
p(t). We refer to [50l [65] for a detailed explanation of the quantified equivalence relations. 
The propagation of chaos for the Vlasov-Poisson equations with singular force has recently 
been investigated in [49]. Here, we are only able to provide such a result in a more restrictive 
setting that in the previous section. Namely, we only show the propagation of chaos for 
d > 3 and with a more restrictive condition on the allowed singularities a > depending 
on the regularity of the initial data 1 < p < oo. 

Theorem 5.1. Given pit) G L°°(0,T; (L 1 nL p )(M d ))nC([0,T],Pi(R d )) the unique solution 

to ([3ZED with initial data p° G Vi(M. d )nL p (R d ), 1 < p < oo, up to time T > 0. Assume that 

p° has compact support, that the initial positions X N >° := {Xf}fL x are iid with law p°, and 

that 

X) — 1 
{l + a)p < - _ - d, 

with a > 0. Then the propagation of chaos holds in the sense that 

P[ sup di(M*),p(i))>-^] -+0, as iV^+oo, 
\te[o,T\ N-y/ d l 
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where 7 is a positive constant satisfying 

p'(2p -!)(! + a) 
dip - 1) 



< 7 < 1. 



Remark 5.1. The condition on a gets more and more restrictive as p gets smaller and 
smaller. In d = 2, even for p = 00 the condition is empty for a > 0. In d = 3, you get 
the condition a < 1/2 for p = 00 and with p = 5 ^ 13 the condition is already empty. We 
also notice that the existence and uniqueness of the solutions are guaranteed by Theorem \4- 1\ 
and taking expectations in the corresponding inequalities for the particle system. Finally, 
in case — 1 < a < 0, the propagation of chaos holds using the same strategy as in Corollary 
\3.1\ by taking expectations in the inequalities for the evolution of the Wasserstein distance. 

We will follow the strategy recently introduced in [19] for the Vlasov equation. We first 
find a deterministic version of the propagation of chaos. This means that we consider a 
regularized system of particles as a kind of middle ground between the solution of the mean- 
field equation (|3,ip and the random particle evolution. More precisely, we define the "blob" 
initial data p^ as 

where e > to be chosen as a function of the number of particles N and Cd is the volume 
of the unit ball in dimension d. We also define the "blob" approximation p^if) to be the 
solution of the system (|3.ip with the kernel W satisfying fl3.5H given by Theorem 14.11 and 
"blob" initial data p® N . 

In the rest, e is chosen as a function of N as e{N) = TV -7 ' with < 7 < 1. It is easy 
to check that ||p^|| p — N^~ l >' p for N large enough, then we can wonder how far is the 
empirical measure to its blob approximation if we assume a bound on ||p^-|| p independent 
of N. 

Proposition 5.1. Under the assumptions of Theorem I5.il and assuming that there exists 
C\ > independent of the number of particles N such that 

\\p%\\p < Ci, and rf m >—e r , 

with 1 < r < i(i +a ) ■ Then, there exists T > such that the solutions pAr(t) and the 
empirical measure pN^t) are well-defined for all t £ [0,T], and 

doc(PN(t),p N (t)) < d 0O (p^,^)e C2T < e(N)e c * T , 
where C<i > is independent of N. 

Proof. We follow a similar argument to Theorem 13.11 We first notice from Theorem 14.11 
that there exists a common time of existence T > of the solutions pn independent of N 
since it only depends on [|pjy|L and a. The empirical measure also exists up to this time 
since it will be smaller than the possible first collision time of particles. Moreover, due to 
(|4.1|) . we get that ||/9jv(t)|| p < C, for all t G [0,T], where C is independent of N. We next 
substitute PN(t) for p(t) in the proof of Theorem 13.11 and thus all estimates in Step A and 
B hold to deduce 



dt 



< c m \\ PN \\ (1 + vf'v m {1+a) ) < c m (1 + 4 P '^ 1+a) ) 
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and 



drj., 
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d/p' -(1+a) 



d /p'~(l+a) 



-C Vm \\ PN \\ [I + r,%* f]^ l+a >) > -C Vm (l + 4 P 'r, m 

where r)^{t) := doo(pjv(*)> //jv(£)). Note that the condition r > 1 makes sense since e ~ 
Vn — ^m — C^ r f° r e small enough. We finally conclude the desired result using a similar 
argument as in Step C of the proof of Theorem 13.11 since 



(ri%) d/p ' 



< Ce d/p '- r{ - l+a) -> as N 



oo, 



□ 



(V?n) 1+a 

by assumption. 

We now present two propositions showing that the assumptions on p® N and 77^ in Propo- 
sition 15.11 are generic in a probability sense when the initial positions X ,0 are iid with law 

2p-l 

p° in L p . We first prove in Proposition 15.21 that r^ n is roughly larger than N d (p-!) if the 
X '° are iid with law p°. 

Proposition 5.2. lei p° G ^(R^) n L p (R d ), 1 <p< 00, and the initial positions X N '° be 
iid with law p°. Suppose there exists L > such that 



2cf \\p°\\ p L7 <N. 



then rffn satisfies 



Proof. Choose an r G 



2p-l 



x° k e 



rf m > LN~^T)\ > e~ 2c d ^°^ LP . 

Then rf^ > r holds if 

l d \A k , with A k = }J B(Xf,r) 



KKfc-l 



2p-l 



for all 1 < A; < TV. It implies from our assumption with r = LN ^p- 1 ) that 

AT 

(x) dx 



rt > LN ^)) > [I 1 - / p% 

k=1 l JA k 



> 



JV-1 

n 

fc=i 






and thus using that ln(l — x) > — 2x if x G [0, ?>], we conclude 



AT-l 



InPfaJJ, > r) > -2c! 7 \\p°\\ p L7 N~ 2+ p J^ k ^ ^ ~ 2c d' \\p°\\p l ^ 



k=i 



□ 



The next proposition gives some bound on the large deviation of ||pjvllp- It states roughly 
that H/o'/vllp is of the same order that ||/o || p , if the X N, ° are iid with law p°. 



Proposition 5.3. Let p° G Vi(M. d ) n L p (M. d ), 1 < p < 00, with compactly support included 
in [— R, R] d . For any iid X N '° with law p°, the smoothed empirical measures p Q N defined 
in (|5.ip satisfy the explicit "large deviations" bound 

HL d \\p% < \\p%\\v) < [2(R + l)] d We- e *W> Nl - r , 



THE DERIVATION OF SWARMING MODELS 27 

where L^ and cr are explicitly given by 

2 In 2 , r 4(4[[v^]] + l) d/p 
cr := j and L d := , 

[2(R + 1)]p Cd 

with [[•]] denoting the integer part. 

Proof. For any Xi S R and x G R , we have 

1 " 1 

8=1 

where ^ stands for the cardinal (of a finite set). Next, we cover [—R,R] by M disjoint 
cubes Cfc of size e d , centered at the points (ck)k<M- The number M of square needed 
depends on N via e, and is bounded by 

d 

M < 



2(12 + 1) 



Assume that x S Cfc for some 1 < A; < M, i.e., |x — Ck\ < -^p, then 

#{? s.t. |x — Xj| < e} < #{z s.t. |cfc — Xj|oo < 2vcfe}, 

and for any 1 < p < oo we obtain 

r d(i-p) 



/ ( P %{x)fdx < ^— 4- #{i s-t. |<* - Xiloo < 2Vdey 
Jc k i N Cd) P 

#{is.t. xGctr, 



Ck 

£ d(l-p) 



(Nc d y 

where C d denotes the cube of center Ck and size (A\/de) d . Let us consider the set of cubes 
of the lattice that contains Ci, i.e., 

Ct C |J C 3 

3 64 

where Ik = {j such that C^ n C-,- ^ 0}. It is direct to check that #1& < Md with M^ = 
(4[[ya]] + l) rf . Moreover, there are only Mj possible values of 1 < k < M such that j £ /& 
for a given 1 < j < M. This yields 

r F d(l-p) M 

/ ^(x)rdx < j— - e E #<* s -t- s G c ^ p 

^-^rS^ 8 -*-^ 667 ^- (5 - 2) 

Let us introduce the notation iVJt := #{i s.t. Xi G C^}. iV& is a random variable which 
follows a binomial law B(N,Sk) with s^ := J c p°(x)dx. If L||/0°||p < ||/o^-|| p , then (|5.2p 
together with Holder's inequality imply that 

M ( At\P M 

E< ^ { -^ L ^\\P°N\\ P P > N?L?e d 7\\ P X > N^L^sl 
fc=i d fc=i 
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where L := c d L/(M dr ) 1 ' p . But, if this happens, it means that for at least one k < M, 

i 

N k > (hr^Nlye^WpX + \ NP ~ LPs i) v 

> -M-pNe7L\\p°\\ p + ^NLs k > — (c R e d \\p% + s fe ) , 

- / 

with cr := 1/ [2(R + 1)] p , where the concavity of x ' p was used. Then, we deduce that 

M ~ 

V(L\\p% < \\p%\\ p ) < J>(iV fe > —[cRe d \\p°\\ p + s k ]) . 

fe=i 

Since N k is a random variable which follows a binomial law B(N,s k ), then for any A, the 
exponential moments of N k are bounded by 

r i N x 

E(e XN ")< l + (e x -l)s k < e (« A -i)^. 
This together with Chebyshev's inequality implies that 

M 

¥{L\\pX<\\p° N \\ p )<Y,ne XNk )e 



M 
.0 . . < n On \ < V- F / \N k ,-\"±[c R e d \\p°\\ p +s k } 



fe=l 

M 

fe=l 



Taking A = InZ/ with the notation L' = ^, we get 



A/ 

-(i'lni'+l-i')iVs fe -L'lni'c i jiVE d ||/9 ||p 



(£||p% < HpaHIp) < 5>" 

fe=l 

M 
< y^ e -L'lnL' Cfl V £ d ||pO|| p = Me ~L> In L'c R Ne d \\p \\ p ^ 



fc=l 

where we used xlnx — x + l>0, for x > 0. With the scaling e(N) = N~d : we get 
V>{L\\p°\\ p < \\p%\\ p ) < [2(R + l)] d N^e^ L ' lnL V^ Nl -\ 

In particular, choosing L = L d = 4(iW^)p jc d so that L' = 2, we get the desired result 

F(L d \\p% < \\pUp) < [2(R + l)] d N^e- c ^°^ Nl ^ , 
for 1 < p < oo. In the case of p = oo, we first notice that as in (|5.2|) . we deduce 



\\Pn\\oo < T7 j sup #{i s.t. c fe - Ailoo < ej = — ? sup N k . 

Nc d e a i<k<M Nc d e a i< k < M 

Since N k follows a binomial law B(N,s k ) and s k < ||p || oo e d , above estimates allow us to 
conclude the desired inequality. D 

We are now in a position to give the proof of propagation of chaos. 
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Proof of Theorem \5.1\ We introduce several sets for the random initial data: 

, , r vNfi . f) \ _r\ . , r yN,0 . r II „0|| v, m .0 || l 

wi .— {A ■r] m >e j, w 2 •— {A . L d \\p \\ p > \\p N \\pi, 
and 

where r, e and L^ are given in Propositions 15. 1 [ 15. 2\ and 15.31 We first provide the estimate 
of ¥(uji). Note that since the assumption on 7, we obtain 

2p-l d 

< 



7(p — 1) p'(l + a)' 

This yields the existence of r verifying 

2 2p - 1 d 

K - < —=- — <r < 



7 j(p - 1) p'(l + a) 

This again implies the existence of /3 > satisfying 

d a 2 P~ 1 

7 lip -I) 

From Proposition 15.21 if we choose L = N~@ , e = N~" t ' d , then 

(wfl = P (A-^ : r& < e r ) = P (X^ : r& < N~^\ 

< P ( X N >° : ^ < LN~^l) < 1 - e- 2c i /P ^°H^ d/p ' 
< 



'D 



2c 1 J p '\\p°\\ p L d /p' <CN~ S , 

for a sufficiently large iV such that N > (2c d \\p°\\p) p ' +dfi , where s = -£■. For the estimate 
of ^((^2)7 we use the result of Proposition 15.31 to obtain 

P(w|) < CiV 7 e- CJv1 " 7 . 
Finally the estimate of P(w|) follows from |15t Proposition 1.2 of Annexe A] (see also 

muni) that 

F(X N >°:d 1 (p N ,p°)>s)<CN~ s ', 

where C and s' are positive constants. We now denote uj := uj\ n W2 H W3. Then we have 

P(w c ) < CN~\ 

for some positive constants C and L If the initial data belongs to w, then we obtain from 
Proposition 15.11 that 

Ce CT 
di(pN(t),PN(t)) < doo(p N (t),p N (t)) < d , for t £ [0,T]. 

We also notice from Theorem 14.11 that 

„0 A „CT s 1 1 / ,,0\ 1 j /,,0 „0An„CT 



di ( P {t),p N (t)) < di (p ,^) e CT < (di (p°,/i^) +d 00 (/&,,&)) 



e 



for all £ £ [0,T]. Since doo ^%,p%) < £ and the initial data belongs to w, this yields 

Ce CT 



di(p(t),pAr(t))< 



N"'/ d ' 
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for all te[0,T] since 

„0 A s- A ( Ji ,.0 \ i a I ,,0 J) 



< 



Ce CT 
N"// d ' 
Hence, we have 



dl {p°,P N ,e) < dl (p°, A) + doc (p N ,P N ,e) 



( Ce CT \ 

P(W) < P [sup dl (p(t), PN (t)) < -j^jj J 



and it implies the desired result 



Ce CT \ C 

™% d iW),PN{t))> w \<nun< w - 
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